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We revisit the power spectrum analysis of the complete sample of the two 
degree field (2dF) QSO redshift (2QZ) survey, as a complementary test of the 
^ work by Outram et al. (2003). A power spectrum consistent with that of the 2QZ 

group is obtained. Differently from their approach, fitting of the power spectrum 
is investigated incorporating the nonlinear effects, the geometric distortion and 
the light-cone effect. It is shown that the QSO power spectrum is consistent 
with the A cold dark matter (CDM) model with the matter density parameter 
Q m = 0.2 ~ 0.5. Our constraint on the density parameter is rather weaker than 
O ■ that of the 2QZ group. We also show that the constraint slightly depends on the 

equation of state parameter w of the dark energy. The constraint on w from the 
QSO power spectrum is demonstrated, though it is not very tight. 
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1. Introduction 



The 2dF QSO redshift (2QZ) survey has established that the QSO sample is a useful 
probe of cosmological models as a tracer of the large scale distribution of mass (Croom, 
et al. 2001;, Hoyle, et al. 2002). In general, constraints on cosmological parameters from 
QSO sample are not very tight. However, the cosmological parameters estimated from the 
QSO sample have a unique implication for cosmology (Outram, et al. 2003; Yamamoto 
2003a). For example, the cosmological principle can be tested by comparing with the result 
from other observations such as galaxy redshift survey and cosmic microwave background 
anisotropies. The pioneering work on the QSO power spectrum analysis was done by Hoyle 
et al. (2002) with the 2QZ 10000 catalogue. Recently Outram et al. have reported the 
result of the similar analysis with the final 2QZ catalogue containing 22652 QSOs (2003). 
They have shown that the QSO power spectrum is consistent with the Hubble Volume 
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ACDM simulation. Furthermore, by fitting the power spectrum with the ACDM model 
within linear theory of density perturbations, they obtained a constraint on the cosmological 
density parameters. 

In the modeling of the QSO power spectrum in Outram et al (2003), however, the 
light-cone effect (Matarrese et al. 1997; Matsubara, Suto & Szapdi 1997; Yamamoto & Suto 
1999), the geometric distortion (Ballinger, Peacock & Heavens 1996; Matsubara & Suto 
1996) and the nonlinear effects (Mo, Jing & Borner 1997; Magira, Jing & Suto 2000) are not 
taken into account. The neglect of these effects might fail to estimate the correct cosmological 
parameters. To test this point, we revisit the power spectrum analysis of the 2QZ sample. We 
have independently performed the power spectrum analysis of clustering with the complete 
2QZ sample. Then we fit the 2QZ power spectrum with theoretical template incorporating 
the effects, which are not considered in the work by Outram et al. (2003). The methodology 
in the present paper is almost same as that in the reference (Yamamoto 2002), in which the 
fitting of the 2dF QSO power spectrum from the 10000 catalogue was investigated using an 
analytic approach. Thus the primary purpose of the present paper is to test the robustness 
of the result by Outram et al. (2003) for independent determination of the power spectrum 
and for more careful modeling of the theoretical power spectrum, including the nonlinear 
effects, the geometric distortion and the light-cone effect. 

On the other hand, Calvao et al. (2002) claimed that the equation of state of the dark 
energy w might be constrained from the 2dF QSO sample. Due to the geometric distortion 
effect, the QSO redshift-space power spectrum may depend on w even if the original matter 
power spectrum (or the transfer function) does not depend on w (Yamamoto 2003b). The 
strategy in the present paper is not the one proposed by Calvao et al. (2002), however, we 
check a constraint on w by considering how the estimated density parameters depends on w 
by the fitting of the power spectrum. The second purpose of this paper is to test the equation 
of state of the dark energy w using the QSO power spectrum. This paper is organized as 
follows: In section 2, we describe our power spectrum analysis. In section 3, our theoretical 
modeling of the QSO power spectrum is explained. In section 4, constraint on the density 
parameters is discussed by fitting the 2QZ power spectrum. Section 5 is devoted to summary 
and conclusions. Throughout this paper we use the unit in which the light velocity equals 1. 



2. The power spectrum analysis 

In our power spectrum analysis, we use the complete sample of the full 2QZ survey, which 
is publicly available (http : //www. 2dfquasar.org/Spec_Cat/). The 2QZ survey covers two 
area of 5 x 75 deg 2 , one in the South Galactic Cap (SGC) and the other in the North Galactic 
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Cap (NGC), respectively, in the range of redshift less than 3. The survey area is defined by 
the equatorial coordinates from a = 21 h 40 to a — 3 h 15 and —32.5° < 5 < —27.5° in the 
SGC, and 9 h 50 < a < 14 h 50 and -2.5° < 5 < 2.5° in the NGC, respectively. The survey 
area of the NGC is jagged and we select a simple rectangle area in our power spectrum 
analysis. Then we use 10713 and 8443 QSOs in the SGC and the NGC, respectively, in the 
range of redshift 0.2 < z < 2.2 1 , incorporating the hole information publicly available. 

We describe the estimator of the power spectrum adopted here. Three dimensional map 
is constructed by introducing the distance 

r dz' 

s(z)= / =, (1) 

Jo iW0-3(l + z'f + 0.7 

where s is the comoving distance of the A cold dark matter model with the density parameter 
0.3. We denote the density field by n(s) and the mean number density by n(s), where s = s^y 
with 7 specifying the direction. Introducing a random synthetic density field n s (s), which 
has mean number density 1/a times that of n(s), we define the Fourier coefficient 

JF(k) = A' 1 ' 2 J ds[n(s) - an s (s)]e !ts (2) 

with A = f cfen(s) 2 . The estimator of the power spectrum is defined 

P(k) = ±rf rfk|f(k)| 2 , (3) 



where 14 is the volume of a thin shell in the k-space with the radius k. In the case a C 1, 
the variance of the power spectrum is 

APW 2 ^. (4) 

Note that we have not used the optimal weighting scheme by setting the optimal weight 
factor being constant (Feldman, Kaiser and Peacock 1994, Tegmark et al. 1998, Yamamoto 
2003b). This choice does not alter the result of the QSO power spectrum analysis because 
the QSO is sparse and nP(k) < 1. Instead of equation (2), the discrete density field can be 
rephrased as 

.F(k) = A-W^ni^e*** - a ^ n s ( Sj )e^], (5) 



1 Thc 2QZ group used the QSOs in the range of the redshift 0.3 < z < 2.2, which is slightly different from 
our choice. This difference does not alter our result. 
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where s, and Sj are the position of the i-th QSOs and the j-th random objects. 

Figure 1 plots the power spectrum (filled squares), which is obtained by combining 
results from the SGC and NGC data sets. Fig.l shows a good agreement with the result 
(open squares) by Outram et al. (2003). However, the error bar of our power spectrum is 
larger than that of 2QZ group. This can originate from the difference of the error estimator 
adopted, which we describe in more details below. The solid curve in Fig.l is the theoretical 
curve of the A cold dark matter (CDM) model with the cosmological parameter, Q m = 0.28, 
tt b = 0.045, h = 0.7, cr 8 = 0.9 and n = 1, motivated by the WMAP result (Spergel et al. 
2003, see also next section). 

Finally in this section, we compare the error estimator of the power spectrum adopted 
in the present paper with that of the 2QZ group. There are two differences. First, they 
used the equation 2.3.2 in paper by Feldman et al. (1994). In the limit that n is a constant, 
there is the discrepancy of the factor 2, if 14, the volume of a thin shell in the k-space, is 
same. In our analysis we used Vk = 4irk 2 Ak, where AA; is the width of the k bin. Our 
error estimator is consistent with that formulated by Tegmark et al. (1998; see equation 
62 in their paper). Second, the 2QZ group adopted the different method for V k) estimated 
by Vk = N k (Ak) 3 with N k the number of independent modes in the /c-shell and (A/c) 3 the 
volume of one /c-mode (for details see also Hoyle, et al. 2002; Hoyle 2000). 



The author considered a constraint on the density parameters, Q m and f^, by fitting 
the 2dF QSO power spectrum found by Hoyle et al. (2002), in a previous investigation 
(Yamamoto 2002). We here adopt the similar methodology to apply it to the 2QZ power 
spectrum obtained in the previous section. The important improvement of the present 
work is the inclusion of nonlinear effects in modeling the theoretical power spectrum and a 
systematic uncertainty in measuring the redshift in the 2QZ survey, as well as the inclusion 
of an arbitrary equation of state parameter w of the dark energy. In this section we briefly 
explain the modeling, following the same notation in the previous paper. 

We restrict ourselves to a spatially flat FRW universe and follow the quintessential cold 
dark matter (QCDM) model. The effective equation of state of the dark energy, w, can 
be a function of redshift, however, we consider the QCDM model with a constant equation 
of state for simplicity (Wang et al. 2000). In this case the relation between the comoving 
distance and the redshift is 



3. Modeling the theoretical power spectrum 
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where H = 100/ikm/s/Mpc is the Hubble constant. Wang & Steinhardt (1998) have given a 
useful approximate formula for the linear growth index in the QCDM model, which we adopt 
in the present work. For nonlinear modeling of the mass density perturbation, we adopt the 
simple fitting formulae for the nonlinear mass perturbation power spectrum presented by 
Ma et al. (1999) for the QCDM model, combined with the transfer function by Eisenstein 
& Hu (1998), which is robust even when the baryon fraction is large. 

To incorporate the nonlinear (Finger-of-God) effect due to the random motion of objects, 
we consider the power spectrum multiplied by the damping factor D^iiOp], where <7p is the 
pairwise velocity dispersion and q\\ is the comoving wave number parallel to the line-of-sight 
of an observer in real space. We adopt the expression assuming an exponential distribution 
function for the pairwise velocity, and the corresponding damping factor is (Mo, et al. 1997; 
Magira, et al. 2000) 

(T) 



The pairwise velocity dispersion is the function of redshift z. We computed a\ by the ap- 
proximate formula for the mean squire velocity dispersion at the large separation determined 
through the cosmic energy equation (Mo, et al. 1997; Magira, et al. 2000; Suto et al. 2000). 

On the other hand, as described in the paper by Outram et al. (2003), systematic 
uncertainty in measurement of the redshift causes an apparent velocity dispersion. We 
incorporate the uncertainty in redshift-measurement Sz, based on the following consideration. 
We write the density fluctuation field <5(s) in a Fourier expansion form as 



S(s) = V5ke !k ' (s+As) , (8) 



D ik-(s+As) 



where As represents the error in measuring the position, which is related with 5z by |As|(= 
As) = 5z(ds/dz). Then an ensemble average of 5(s)5(s') can be written as 

(5(s)5(s>)) = ]T(|<5 k | V Ms - s,) (e* k -( As - As ')) , (9) 

k 

where we have assumed that 5k and As are independent probability variables and (5 k ^k') = 
(|5k| 2 )5® (k — k'). We further assume that the angular coordinates are well determined and 
As = 7 As. In this case, denoting the wave number of the line-of-sight direction by k\\, we 
have 

k.(As-As')\ = e -^((A s -As') 2 /2) = e -^(As 2 ) ; ^ 
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where we have assumed that As and As' are independent Gaussian probability variables with 
the variance (As 2 ). Therefore the damping factor due to the error in redshift-measurement 
is written 



where (Sz 2 ) is the variance of the error in measuring the redshift. In the present paper 
we adopt 5z = 0.0027,2 following Croom et al. (2003). Figure 2 shows the damping fac- 
tor _D[/c||<xp] (solid curves) in the ACDM model 2 , as function of the wave number k\\ with 
the redshift fixed as z = 0.5, z — 1 and z = 1.5 for three curves from top to bottom. 



considering the angular averaged power spectrum, and this figure does not exactly express 
the damping factor of our power spectrum. However, this figure indicates that the damp- 
ing factor have a substantial effect on the power spectrum shape. Thus the uncertainty in 
redshift-measurement can be an important factor in the power spectrum analysis. 

Concerning the modeling of the bias, we assume the scale independent model. Following 
the previous work (Yamamoto 2002), we consider the model b(z) = b /Di(z) v , where b and 
v are the constant parameters and D\(z) is the linear growth rate normalized as -Di(O) = 1. 
We determine bo to minimize x 2 , which we define in the next section. In the present paper, 
we show the result adopting v — 1. Our result slightly depends on u, however, plausible 
alternation of v does not alter our conclusion qualitatively, (see also Yamamoto 2002). We 
consistently use the number density per unit redshift per unit solid angle dN/dz from the 
catalogue. 

It will be useful to show how the theoretical power spectrum depends on the equation of 
state w of the dark energy. Curves in Figure 3 are theoretical power spectra with the various 
equation of state w = —1 (solid curve), w = —2/3 (long dashed curve) and w = —1/3 (dashed 
curve). The other cosmological parameter is fixed as VL m = 0.28, Vl b = 0.1, h = 0.7, er 8 = 0.9 
and n — 1. The difference of the three curves comes from the geometric distortion effect: 
The comoving distance in real space is given by equation (6), while the power spectrum is 
measured in redshift space with the radial coordinate s(z) defined by equation (1). Therefore 
the power spectrum in redshift space is distorted compared with that in real space. Denoting 
the wave number in real space by q, the geometric distortion is described by the scaling of 
wave numbers q k/c, where c is the ratio dr/ds and r/s for the component parallel and 
perpendicular to the line-of-sight direction, respectively (for details see e.g., Ballinger et al. 




(11) 



The dashed curves show the damping factor D[k\\ap] multiplied by X>[fo;]. Because we are 



2 For definiteness, D[k\\<j-p] denotes the damping factor D[<7||crp] with replacing g|| with (ds/dr)k\\. 
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1996, Yamamoto 2002; 2003b). For the model with w > — 1, the factor c becomes less than 
unity. This shifts the power spectrum from right (large k) to left (small k) for the model 
with w > — 1, as is clearly shown in Fig. 3. Thus w can be measured from the redshift 
power spectrum due to the scaling effect of the geometric distortion. Such a test using the 
scaling effect is distinguished from the Alcock-Paczynski test which relies on the anisotropy 
of clustering in redshift space (e.g., Ryden 1995; Alcock & Paczynski 1979). 



4. Discussions 

In this section we compare our observational power spectrum from the 2QZ catalogue 
with the theoretical model to constrain cosmological parameters. We introduce 

_ " [P^fcQ - P^f 
x AP(^) 2 ' { ' 

where P th (ki) is the value of the power spectrum from the theoretical model at the wave 
number ki, P oh {ki) is the observational data and AP(fcj) is the variance of errors, for which 
we use the (N =)20 data points obtained in Fig. 1. Following the analysis by Outram et al. 
(2003), we consider the likelihood function oc e~ x2//2 . Figure 4 shows contours of the likelihood 
function with \ 2 found when different values of the cosmological density parameters, Q m and 
Q{,/Q m are used in the theoretical modeling. Each panel corresponds to constraints for the 
difference model of P th (k): The panels (a), (c) and (d) assume w = —1, —2/3 and —1/3, 
respectively. The panel (b) is same as (a), but switched off the damping factor due to the 
error in the redshift-measurement by setting Sz = in equation (11). In Fig. 4, we fix the 
other parameters h = 0.7, n = 1 and <r 8 = 0.9. The curves are contours of confidence 65%, 
95% and 99% on the plane. 

First, the panel (b) in Fig. 4 shows the effect of the error in the redshift-measurement in 
modeling the power spectrum. The peak value of (a) is located at Q m = 0.27, while that of 
(b) is located at Q m = 0.24. Thus the effect of 5z can be of importance for determining fl m , 
precisely. Second, Fig. 4 demonstrates the constraint on the density parameters depends on 
the equation of state parameter w of the dark energy. It shows that the preferable value 
of Q m and Q b increases as w becomes large. Figure 5 shows the contour of the likelihood 
function with x 2 found when different values of Q m and w are used in the modeling the power 
spectrum, where we fixed f2 b = 0.045, motivated from the WMAP result (Spergel et al. 2003). 
It is shown that the density parameters in the range 0.2 £ Q m £ 0.4 is preferable, which is 
not very sensitive to the value of w. Therefore the constraint on w is not very tight, and 
w £ — 0.2 (—0.1) is only excluded at the one (two) sigma level. 
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5. Conclusion 

In summary we have revisited the power spectrum analysis of the complete 2dF QSO 
sample, as a complementary test of the work by Outram et al. (2003). Our analysis has 
reproduced a power spectrum, which is consistent with the 2QZ group. We have investigated 
the fitting of the power spectrum including the light-cone effect, the geometric distortion and 
the nonlinear effects. It is shown that the QSO spatial power spectrum is consistent with 
the ACDM model with Vl m = 0.2 ~ 0.5, which is effectively consistent with the 2QZ group. 
However, our constraint is weaker than that of the 2QZ group. This will be traced back to 
the difference of the error estimation of the power spectrum. In the present paper, we have 
emphasized the importance of the error in measurement of the redshift 5z because it may have 
influence in estimating the density parameters. We have also investigated the effect of the 
equation of state of the dark energy w. For models with large value of w > —1, the preferable 
value of the density parameters becomes large. When we fix the baryon density parameter as 
that found by the WMAP team, the matter density parameter 0.2 <, VL rn <, 0.4 is preferable. 
No tight constraint on w is obtained, w >, — 0.2 (—0.1) is only excluded at the one (two) 
sigma level. In general QSO sample is sparse and the shot noise contribution is substantial. 
Therefore the constraint from it is not very tight. However, such a power spectrum analysis 
will provide more useful constraint on w when applied to the SDSS luminous red galaxy 
sample (Yamamoto 2003b) and next generation redshift survey such as the KAOS project 
(Seo and Eisenstein 2003). 

The author thanks anonymous referee for useful comments on earlier version of the 
manuscript, which helped improve it. This work is supported in part by Grant-in- Aid for 
Scientific research of Japanese Ministry of Education, Culture, Sports, Science and Technol- 
ogy, No. 15740155. 
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Fig. 1. — Power spectrum from the 2dF QSO sample. The filled square is the result of 
our analysis, while the open square is from Outram et al. (2003). The solid curve is the 
theoretical curve of the ACDM model with the cosmological parameter, Vl m = 0.28, VL b = 
0.045, h = 0.7, a$ = 0.9 and n — 1. (see also section 3) 
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Fig. 2. — Damping factor D[k\\ap(z)] (solid curves) as function of the wave number k\\ with 
the redshift fixed z = 0.5, z = 1.0 and z = 1.5 from top to bottom. The dashed curves show 
the damping factor multiplied by £>[fc]. We assumed the same theoretical model as that in 
Fig. 1. 
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Fig. 3. — Theoretical power spectra with w — — 1 (solid curve), w = —2/3 (long dashed 
curve) and w = —1/3 (dashed curve). The other cosmo logical parameter is fixed as Q m = 
0.28, Qb = 0.1, h = 0.7, o"8 — 0.9 and n — 1. The amplitude is normalized as P(k) = 
10 4 (/i _1 Mpc) 3 at k — O.l/iMpc -1 . Note that the model of large baryon fraction is adopted 
to emphasize the scaling effect of the geometric distortion. 



-14- 




Fig. 4. — Contours of the likelihood function of our QSO power spectrum. The panel (a), 
(c) and (d) assumes w = —1, —2/3 and —1/3, respectively, in which we adopt h = 0.7, 
n = 1 and og = 0.9. The panel (b) is same as (a) but neglected the damping factor due to 
the error in redshift-measurement. In each panel, contours are confidence of 65%, 95% and 
99%. The cross point shows the best fitted parameter: (Q m ,Qi,) = (0.27,0.06), (0.24,0.06), 
(0.33,0.08), (0.41,0.11), from (a) to (d), respectively. The dashed curve shows Q b = 0.045 
for comparison. 
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Fig. 5. — Contours of the likelihood function on the Vt m and w plane. Contours are confidence 
of 65%, 95% and 99%. Here we assume tt b = 0.045, h = 0.7, n = 1 and a 8 = 0.9. The 
dashed line is Q m = 0.28. 



